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Abstract 

Multi-walled hollow nanoparticles made from tungsten disulphide 
(WS2) show exceptional tribological performance as additives to liq- 
uid lubricants due to effective transfer of low shear strength material 
onto the sliding surfaces. Using a scaling approach based on contin- 
uum elasticity theory for shells and pairwise summation of van der 
Waals interactions, we show that van der Waals interactions cause 
strong adhesion to the substrate which favors release of delaminated 
layers onto the surfaces. For large and thin nanoparticles, van der 
Waals adhesion can cause considerable deformation and subsequent 
delamination. For the thick WS2 nanoparticles, deformation due to 
van der Waals interactions remains small and the main mechanism 
for delamination is pressure which in fact leads to collapse beyond a 
critical value. We also discuss the effect of shear flow on deformation 
and rolling on the substrate. 

1 Introduction 

Graphite and layered material made from metal disulphides (M0S2, WS2) 
and similar composites (MoSe2, WSe2, BN, etc.) are good lubricants since 
under shear the layers can easily slide over each other due to atomic smooth- 
ness and weak van der Waals (vdW) interactions JIJ. Therefore they are 
widely used as solid lubricants or additives to liquid lubricants. However, 
finite sized crystallites in powders have edges with dangling bonds which can 
react chemically with the sliding surfaces. This problem can be avoided by 
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using hollow nanoparticles made from the same material, where the layers 
are not planar, but are bent into closed shells. It has been shown recently 
that multi- walled WS2 nanoparticles perform very well as additives 0, I3J. 
Closer investigations using the Surface Force Apparatus (SEA) revealed that 
WS2 layers are transfered onto the sliding surfaces by delamination of the 
nanoparticles 0]. These peeled layers form islands on which the Friction 
Force Microscope measured much smaller friction than on the surrounding 
mica. Hence multi-walled nanoparticles can act as reservoirs which release 
layers of low shear strength exactly where needed while avoiding too many 
dangling bonds. 

A theoretical understanding of the mechanical properties of multi-walled 
hollow nanoparticles is important for future tribological applications. In 
this paper, we use scaling arguments based on continuum elasticity theory 
and pairwise summation of vdW interactions in order to investigate theo- 
retically the effect of adhesion, pressure and shear flow on the deformation 
and mechanical stability of hollow nanoparticles of spherical shape. Figure 
[1] schematically depicts some of the aspects discussed in the following. We 
show that vdW adhesion to the substrate (as well as to each other) can be 
several orders of magnitude larger than thermal energies and scales linearly 
with radius, but is almost independent of thickness. Despite the large en- 
ergy of adhesion, we find that for typical WS 2 nanoparticles, coherent and 
even incoherent deformations due to the vdW adhesion to the substrate are 
only in the Angstrom range and will leave the nanoparticles basically intact. 
Thus for WS2 nanoparticles the vdW adhesion favors the release of delami- 
nated layers onto the surfaces, but does not trigger the delamination itself. 
The main mechanism for delamination of WS2 nanoparticles is shown to be 
pressure which leads to a mechanical instability in linear elasticity theory. 
This might explain why damage of WS2 nanoparticles was found to occur 
only beyond a critical load ||. We also show that shear flow does not lead 
to considerable deformation of WS2 nanoparticles and that very large shear 
rates are needed in order to make adhering particles roll in shear flow. 

Similar methods have been used before to predict the shape of fullerenes 
0, || and to account for the faceted shape of metal disulphide particles || . 
Compared with ab initio methods ||, tight-binding schemes 0] and molecu- 
lar simulations which have been used before to investigate the mechanical 
properties of fullerene-like material, our approach has the advantage that it 
is asymptotically correct for large systems and universal in the sense that dif- 
ferent material systems enter on the level of their elastic and vdW constants. 
In particular we discuss C, M0S2 and WS2, although we focus on the WS2 
nanoparticles used in the tribological experiments mentioned above. Our 
scaling approach allows us to predict how tribological performance depends 
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on radius R, shell thickness h and layer thickness a of the nanoparticles. We 
show that scaling with h does not result from the vdW contributions, but 
rather from the scaling of the elastic constants with h. Several critical quan- 
tities derived scale strongly with R/h for coherent deformation and with R/a 
for incoherent deformation. Hence tribological performance can be tuned by 
varying the concentration of defects (which switches between the coherent 
and incoherent regime) and the ratios R/h and R/a, respectively. 



2 Preliminaries 

In the following we consider hollow nanoparticles with outer radius R and 
thickness h which results from nesting several elastic shells of thickness a 
each. Typical WS2 particles used in tribological experiments have R ~ 60 
nm and h fa 9.3 nm (15 layers with an interlayer distance of a = 0.62 nm, 
R/h « 6, R/a fa 100). In fig. ^| a high resolution transmission electron 
micrograph of a WS2 nanoparticle with 11 layers is shown. The elastic shells 
are closed since their formation is driven by the energy reduction due to the 
absence of dangling bonds. In order to achieve closure, topology dictates that 
for carbon shells 12 pentagons have to be inserted into the graphite network 
of hexagons. For hexagonally layered composites, the sheets have a more 
complicated molecular structure (typically triple layers) and different kinds 
of defects have to occur to ensure closure. In our continuum approach, we 
assume that defects are distributed in a homogeneous way to give a spherical 
shape with a certain preferred radius R @. Similar approaches to fullerene- 
like particles have considered sheets which are planar in equilibrium; in order 
to account for their curvature, they explicitly considered defects [|, || [7[ §. 
However, in this approach the elastic response of a spherical shell is very 
complicated and no unifying scaling approach is possible. 

The new feature of the elasticity of a closed shell with a preferred radius R 
is that stretching is a first-order effect and the spherical shell cannot be bent 
without being stretched |~2] . This interplay between bending and stretching 



has been studied before for thermal fluctuations of polymerized vesicles |13 
The elastic behavior of the shell is determined by two contributions: bending 
energy Ej, ~ k J dA c 2 with bending rigidity k and mean curvature c, and 
stretching energy E s ~ G J dA e 2 with in-plane stretching modulus (or two- 
dimensional Young modulus) G and in-plane strain e. Here / dA represents 
the surface integral; the elastic contributions in shell theory follow by inte- 
grating over the thickness of the shell. Consider a spherical shell of radius R 
which is expanded by AR. Then c changes by AR/R 2 and e by AR/R, thus 
both contributions are first order effects. Since Ej,/E s ~ k/GR 2 , the rela- 
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tive strength of bending and stretching depends both on material parameters 
and radius; on large length scales R 3> (k/G) 1 ^ 2 and stretching will always 
dominate. 

In the framework of continuum elasticity theory, the elastic moduli k and 
G can be calculated from the in-plane stretching elastic constant Cn of the 
corresponding hexagonal layered material as k = Cnh 3 /12 and G = Cnh 
0. The values for C u are 1060, 238 and 150 x 10 10 erg/cm 3 for C, MoS 2 
and WS2, respectively |TJ]]. The same scaling with h is found for thin films 



made from isotropic elastic material [12 . It only applies for coherent bending 



of the different layers; for incoherent bending, slip occurs between adjacent 
layers and overall bending becomes easier. The bending rigidity then follows 
as k ~ Cna 3 (h/a) = Cna 2 h where a is the effective thickness of a single 
layer . It should be of the same order of magnitude, but somewhat smaller 
than the interlayer distance (which is 3.4 A for C and 6.2 A for M0S2 and 
WS 2 ). Note that the in-plane modulus G ~ Cua(h/a) = Cnh stays the 
same since it scales linearly with thickness. The change from coherent to 
incoherent bending can be considered to be defect mediated; dislocations will 
proliferate for increasing thickness and finally will lead to grain boundaries 
which account for the typical faceted shape seen in fig. M . 



3 Adhesion 

We consider a substrate which interacts by attractive vdW forces with a film 
of thickness h which is a distance D away from the substrate. Pairwise inte- 
gration of the potential — A/-K 2 r 6 (where the Hamaker constant A is typically 
of the order of 10~ 12 erg) over the two volumes yields the vdW energy per 
unit area 

= A h(h + 2D) 
U 12irD 2 (h + DY ' 1 ] 

Here D is an atomic cutoff for the vdW interaction which in the following is 
chosen to be 1.65 A ]Tj|. Although the adhesion energy u scales linearly with 
h for small h, for more than one layer we have h > D and u saturates at a 
constant value u = AjYli^D 2 w 100 erg/cm 2 . For the case of a single layer, 
one can replace the volume integral over the film by a surface integral times 
the effective thickness h; this is equivalent to considering h <C D in eq. (|l|) and 
leads to u — Ah/QirD 3 . Since D and h are of the same order of magnitude, 
we essentially recover the result for a multi-layered film. Hence we conclude 
that due to the rapid saturation of the vdW binding energy with increasing 
h, for multi-layered films under adhesion conditions no considerable scaling 
with h is expected from the vdW terms. 
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It can be shown by similar arguments (e.g. by using eq. (|I|) within the 
Derjaguin approximation) that the energy of a hollow nanoparticle adhering 
to a substrate does not depend on h as long as h > D. Then the adhesion 
energy equals the vdW energy of close approach between a sphere and a 
substrate, Ea = AR/6D |0|. Using values for A, D and R as given above 
yields Ea = 6 x 10 -11 erg = 1400 kT. The energy of adhesion between two 
particles is only a factor of 2 smaller than the one between a particle and the 
substrate, and hence well above thermal energies as well. We thus conclude 
that the vdW interaction leads to considerable adhesion of the particles to 
the substrate and to each other. The same holds true for peeled off outer 
layers: with R = 60 nm, the energy of adhesion can be estimated to be 
R 2 u = 10 5 kT. 



4 Deformation under adhesion 

The onset of delamination can be estimated by considering large deforma- 
tions of the nanoparticles. We do not offer a theory for fracture and estimate 
the limit of mechanical stability by an internal criterion within our contin- 
uum description: fracture and delamination sets in when so much stress has 
accumulated that the deformation becomes of the order of the radius. De- 
termining the deformation of a spherical elastic shell is a difficult problem 
which depends on prestress, elastic constants and size as well as on the nature 
of the deforming force. Here we discuss the deformation on a substrate for 
two cases which can be treated in the framework of our scaling approach. In 
the following H is the indentation. For small deformations H < h (fig. |3]a), 
the shell flattens at the bottom. For large deformations H > h (fig. [|b), a 
contact disc with radius r develops and the elastic energy is localized in the 



circular fold surrounding it. The respective elastic energies are [12 



G l/2 K l/2 ^1/4^3/4 
E sm all ~ „ H , Ei ar g e ~ H . (2) 

ti ri 

Here we neglect a bending term ~ kH/R. We also assumed that due to low 
friction, slip can occur between the contact disc and the substrate; otherwise 
compression energy ~ GH 3 /R would accumulate in the contact disc. The 
rest of the relaxing shell is assumed to keep a spherical form with radius 
R 1 > R in order to keep area constant and to avoid compression energy. 
However, the difference is (R' — R) ~ H 2 /R and can be shown to be negligible 
in the following. 

If the deformations are driven by vdW adhesion, the energy gained on 
adhesion can be shown to scale as Ea ~ uRH for both small and large 
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deformations. Setting this adhesion energy equal to the estimates for the 
elastic energies from eq. (0) yields estimates for indentation: 



R 2 u ( u \ / R \ R A u 2 ( u \ 2 fR\ 3 

H small ~ ~ j j J2, fl^ ~ ~ j j 

(3) 

Here and in the following we will always give the result both in terms of the 
elastic moduli and their scaling with thickness for coherent bending. The 
dimensionless quantity u/Cuh is the ratio of van der Waals adhesion energy 
to two-dimensional Young modulus and will be of order 10 -4 for the WS2 
nanoparticles. In this case we are in the regime H < h and deformations 
will be in the Angstrom range. Thus delamination cannot be expected to be 
caused by adhesion. For incoherent bending, R/h has to be replaced by R/a 
in both regimes. 

One system for which vdW adhesion leads to strong deformations are 
single-walled hollow nanoparticles. In this case, care has to be taken to use 
the correct values for G and k. Since G is a purely two-dimensional quantity, 
G = Cuh can be used where h is identified with the interlayer distance. 
However, k has to be extracted from molecular calculations [|, [10], [IT| . For 
carbon, one finds G = 3.6 x 10 5 erg/cm 2 and k ~ 1.6 x 10~ 12 erg fh 40 
kT. For both M0S2 and WS2, G is smaller by a factor 4 and k is larger by 
a factor 10. In all three cases, for radii larger than a few nm considerable 
deformations arise. The critical radius R c where indentation H and radius R 
become of the same magnitude follows from the second relation of eq. @ as 
For C, we find R c ~ 5 nm and for the metal disulphides 
R c 10 nm. For nanotubes, we find similar results which are in good 
agreement with experiments and molecular calculations [T6|, [T7 . 



5 Deformation under pressure 

We now proceed to show that pressure can lead to an instability which could 
account for the observed delamination. We consider the formation of a con- 
tact disc for a shell which is pressed onto the substrate by pressure in the 
surrounding liquid. The corresponding energy is E p ~ —pRH 2 , thus it scales 
more strongly with H than the restoring elastic energy Ei arge in eq. (0). This 
scaling indicates an instability: small deformations are suppressed, but large 
deformations grow without limits beyond the critical indentation H c where 
Eiarge + E p attains a maximum [[12[ . Setting H c = R gives an estimate for 
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the critical pressure for delamination: 



£1/4^3/4 
R 5 / 2 




(4) 



Using typical values gives p c ~ 1.7 GPa. Incoherent deformations can be 
treated by replacing (h/R) 3 ^ 2 by (a/R) 3 ^ 2 in eq. <Q) and hence would decrease 
the estimate to 25 MPa. Although a full treatment of the instability requires 
a theory for the non-linear elasticity, we can still suggest that pressure will be 
a very likely mechanism to cause delamination. This prediction agrees nicely 
with the experimental observation that damage sets in at a finite threshold 
of the load ||. Note that for single-walled buckyballs C 60 with R = 3.55 A, 
eq. d|]) gives p c = 15 GPa which is in surprisingly good agreement with the 
value 20 GPa found experimentally |TS . 



Although the pressure values estimated here seem to be high, they re- 
fer to particles of nm-size, and the corresponding forces are of the order 
of 10~ 3 dyn. In friction experiments, the pressures and forces needed for 
delamination are most likely to occur near asperities. This might explain 
why delamination has been observed in the SFA only when shearing the two 
surfaces 0]: the nanoparticles then behaves as granular material which gets 
jammed and pressure and forces become localized. 



6 Shear flow 

Since the WS2 nanoparticles are immersed in liquid lubricant in strong mo- 
tion, we finally consider the effect of shear flow. Balancing the deforming 
viscous force with the elastic restoring force and setting the resulting de- 
formation equal to the radius gives an estimate for the shear rate at which 
delamination sets in. For h < R, we find that only stretching is relevant and 
that the critical shear rate is 

G Cu h 
r/R 77 R 

Using viscosity rj w 1 cP (which applies both to water and mineral oil), we 
find 7 C ~ 2.5 x 10 13 Hz. Such high shear rates are unattainable even near 
contacting asperities in macroscopic friction experiments, thus delamination 
is not likely to occur due to shear flow alone. 

Shear flow past nanoparticle adhering to a substrate can cause them to 
roll. The particle begins to roll if the viscous drag evaluated at the midpoint, 
Fs = 6111] R 2 j, equals the friction force Fp caused by the adhesive load 



7 



Fa = 2nRu. Assuming Amontons' law for rolling friction, Fp = [iFa, we 
find 

7c = 3^ (6) 

where the coefficient of rolling friction /i can be taken to have the typical 
value 10 -3 ; then j c is of the order of 5 x 10 5 Hz. This value is higher 
than typical shear rates probed with the SFA (v = 10 /im/s, 7 = 10 4 Hz 
on nm-separation), but certainly lower than shear rates occurring near the 
asperities of macroscopic friction experiments (v = lm/s, 7 = 10 9 Hz on 
nm-separation) . 



7 Conclusion 

We have shown that in friction experiments with multi- walled hollow nanopar- 
ticles, delamination is likely to be caused by pressure. When nanoparticles 
delaminate, dangling bonds develop as they do for powders from layered 
material. However, there are two reasons why they should be less problem- 
atic. First delamination takes place with one or two layers and the density 
of dangling bonds remains low. Second we showed that delamination may 
occur under adhering conditions. The adhesion (and possibly the shear) will 
cause the layers to align parallel to the substrate, in contrast to the powders 
where the orientation is preferentially perpendicular to the substrate. It has 
been shown experimentally that delaminated layers from nanoparticles form 
islands while platelets form a thicker, more disorganized film ||. Our scaling 
approach shows that there are two control parameters which might be used 
for tribological optimization: the ratio of radius to thickness and the con- 
centration of defects (that determine the regimes of coherent and incoherent 
bending). 

Using nanoparticles as additives to liquid lubricants raises interesting 
questions about the role of wetting. The lubricant used in the friction ex- 
periments mentioned is a bad solvent for the nanoparticles, which provides 
an additional force to drive them onto the substrate. A dewetted region 
is expected to form between the adhering particle and the substrate; like a 
capillary bridge in a good wetting situation, it might enhance adhesion. 
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Figure 1: Several aspects of friction experiments with mult i- walled nanopar- 
ticles: (a) weak deformation due to adhesion onto a substrate, (b) strong 
deformation due to pressure, (c) deformation in shear flow and (d) rolling of 
adhering particles in shear flow. 
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Figure 2: High resolution transmission electron micrograph of a multi-walled 
WS2 nanoparticle. The slightly faceted shape is typical and can be attributed 
to the assembly of defects into grain boundaries 0. 
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Figure 3: Possible deformations of a spherical elastic shell on contact with 
a substrate: (a) flattening and (b) contact disc. R is the initial particle 
radius, R! the radius after flattening, h the shell thickness, H the indentation 
and r the radius of the contact disc. For small H, one finds geometrically 
r 2 = 2HR'. 
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